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Abstract 

We calculate the electromagnetic corrections to the decays rj — > 3ir in next to 
leading order in the chiral expansion. We find that the corrections are small in 
accordance with Sutherland's theorem and modify neither rate nor the Dalitz plot 
distributions noticeably 
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1 Introduction 



In the isospin limit m u = ma and e = 0, the decay 77 — > 37r is forbidden. Conse- 
quently, the decay amplitude receives contributions from the QCD isospin violating 
interaction 

n QCD (x) = {md ~ mu) (dd - uu)( x ), (1.1) 

and from the electromagnetic interactions 

^qed(x) = -\e 2 J dyD^(x - y)T(j^(x)j v (y)), (1.2) 

where D^ v {x — y) is the photon propagator. Using soft pion techniques, it has been 
shown that the electromagnetic contribution is much too small to account for the 
experimentally observed rate [1], §]. Much work has been devoted to study the effects 
of the QCD contribution, equation (|1.1|). Gasser and Leutwyler || carried out the 
one-loop calculation within chiral perturbation theory ||, |7)- Although they 
found large corrections to the rate their result failed to reproduce the experimental 
rates |J. Recently, the complete unitary corrections were evaluated 0. Despite 
considerable uncertainties, also these corrections are not quite sufficient to account 
for the observed rate, if the usually assumed value of {m u —m d ) |Ll]] is taken. Instead, 
a somewhat smaller up-quark mass is required. 

In this situation, it is of interest to reconsider the electromagnetic corrections. 
Strictly speaking, they may be divided into the indirect ones, affecting the param- 
eters such as F n which enter the calculation and the direct ones, specific for the 
process under consideration. In the following we will be concerned with the latter 
only. Corrections to F n were calculated previously |TTJ and change its value by about 



1 %, leading to a noteworthy increase of the eta decay rate by about 4 %. 

Within the framework of chiral perturbation theory, electromagnetic correction 
can also be described by a series of effective operators of increasing power in mo- 



mentum or masses of the mesons |12[ . Sutherland's theorem ||] states that the first 
correction of order e 2 p° vanishes where p is a typical momentum. Thus, the electro- 
magnetic corrections are at most of order e 2 p 2 . If the decay amplitude is assumed 
to be linear in the Mandelstam variables, these corrections are further suppressed, 
i.e. p 2 is of order M 2 rather than M 2 0, |13| . 

A similar result (Dashen's theorem [Q) for the electromagnetic mass differences 
states that they are equal for the pions and for the Kaons to leading order in the chi- 
ral expansion. Recently, it was argued that this equality could be violated |l5j (see, 



however (T^|), and one might expect then that also Sutherland's theorem receives 
larger corrections than previously thought which could enhance the decay rate of 
the rj. In this article we determine the 0(e 2 p 2 ) electromagnetic corrections to the 
decay amplitudes 77 — > 37r. We show that they can be safely neglected compared 



3 The Dalitz plot distribution does not get equally large corrections from higher orders in the 
chiral expansion. It is therefore possibly more sensitive to electromagnetic corrections. 
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to the QCD isospin violating contributions, both for the rates and the Dalitz plot 
distribution. The paper is organized as follows: in sections 2 and 3 we review the 
chiral lagrangian for electromagnetic interactions up to order p 2 a. In section 4 we 
calculate the electromagnetic corrections to rj — > 3ir to one-loop chiral perturbation 
theory. The numerical results are given in section 5 and section 6 contains our 
conclusions. 



2 Effective low-energy Lagrangian 

To lowest order in Chiral Pertubation Theory (ChPT), including the electromagnetic 
interactions, the Lagrangian £ 2 is given by ||, [7|, [lTj 

+ ^±(V,U^W) + ^(U^ + X^) + C(Q R UQ L W). (2.3) 

Here, f^ u is the field strength tensor of the photon field 

fnv = dfj,a u - d u a^. (2.4) 

The parameter £ is the gauge fixing parameter, which is set to £ = 1 henceforth. 
The pseudoscalar meson fields are contained in the usual way in the matrix U, and 
the field x incorporates the coupling of the pseudoscalar mesons to the scalar and 
pseudoscalar currents s and p 

X = 2B (s + ip) =2B M + ... (2.5) 

where the quark mass-matrix is 

/ m v 



M 



m d I . (2.6) 
\ m s 



The covariant derivative V M C/ defines the coupling of the pseudoscalar mesons to 
the photon field a M , the external vector current V 1 , and the external axial current 
A 11 

V"C/ = d"U - i(V + Q R a" + A^)U - iU(V + Q L a" - A"). (2.7) 



The two spurios fields Qr and Ql are introduced to construct a SU(3)l ® SU(3)r 
invariant Lagrangian. In the later calculation they will be fixed by Q R = Ql = Q, 
where Q is the charge matrix of the three light quarks 



(2.8) 
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The constant F is up to order 0{m q ) the pion decay constant. 

The low energy constant C determines the electromagnetic part of the masses 
M n ± , Mk± i of the charged pion and Kaon in the chiral limit 

C( QUQW) = -^(tt+tt- + K+K~) + .... (2.9) 

7T 

Clearly, C( QUQW) contributes equally to the square of the masses of tt^ 1 and K ± , 
in agreement with Dashen's theorem fl4H . 

Since the external currents V 1 and are of order 0(p), where p is their external 
momentum, the product Qa^ must have dimension 0(p). In order to formally 
maintain consistent chiral counting, it is convenient to set [0 

e ~ 0(p), a» ~ 1, (2.10) 

such that C 2 is 0(p 2 ). 



3 Effective Lagrangian to order 0(p A ) 

At order 0(p 4 ) the generating functional has three different types of contributions: 

• An explicit local action of order 0(p 4 ) 

• The one-loop graphs associated with the lowest order Lagrangian £ 2 . 

• The anomaly is of order 0(p 4 ); its contributions will not be discussed here. 
For a discussion, the reader is referred to the literature [0, [7|. 

Since we are not concerned here with the strong infractions,^], we only give the 
electromagnetic terms. The most general chiral and Lorentz invariant, P and C 
symmetric Lagrangian of order 0(p A ) is |12] 

+ K 3 F 2 {{V^QU)(V^QU) + (WUQW)(X7^UQW)) 

+ K 4 F 2 ( V»tfQU){VJJQlfl) + K 5 F 2 { {Vtft, V M f/} Q 2 ) 

+ K 6 F 2 ( VWVpUQWQU + VUV^WQUQW) 

+ K?F 2 ( X ^U + W X ) ( Q 2 ) + K 8 F 2 ( X ] U + W X ) ( QUQW) 

+ K 9 F 2 ( (x ] U + W X ) Q 2 + ( X W + U X ] ) Q 2 ) 



+ K 10 F 2 ( (x ] U + W X ) QU^QU + ( X W + Ux ] ) QUQW) 

+ K U F 2 ( (x ] U - W X ) QU^QU + ( X f/t _ c/ x t) QUQW) 

+ K 12 F 2 ( V»W [V M Q R , Q]U + WU [V M Q L , q] W) 

+ K 13 F 2 (\7»Q R U\7^Q L W) + K 14 F 2 ( V^Q R V ^Q R + WQ L V„Q L ) 

+ K 15 F*( QUQW} 2 + K W F*( QUQW)( Q 2 ) + K 17 F*( Q 2 } 2 , (3.11) 
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where F^ v and F^ v are the field strength tensors of Fff and of F£ respectively 



F L 



V„ + Q% + A p 
Vn + Q L a ll - A. 



F 1 = 8 F 1 — 8 F 1 - 
and the covariant derivative V^Q 1 is defined as 



F 1 F 1 

± 111 ± V 



R,L 



V^Q 1 = d li C? + i\Q I ,F! 



R, L. 



(3.12) 



(3.13) 



Under a chiral transformation the covariant derivative V^Q 1 transforms like the 
corresponding Q 1 



V M Q 7 

V R , L 



SU(3) R , L 



R,L 



(3.14) 



We have set Q R = Q L = Q. Only in the covariant derivative V^Q 1 have we kept 
the index / = R, L in order to make explicit the proper chiral transformation. 
The coupling constants K\ . . . Kn are not determined by chiral symmetry. While in 
principle calculable from the fundamental lagrangian of QCD and QED, we consider 
them as phenomenological constants. K±, Kj, K\q are electromagnetic corrections 
to Fq, B q and C respectively. 

In the next step, the loops generated by the Lagrangian £ 2 are calculated. They 
lead to terms proportional to p A . Their divergencies are cancelled by infinite coun- 
terterms in the constants Ki. The procedure is standard and will not displayed here 
in detail | 12fl . Using the usual dimensional regularization scheme with scale /x, one 
finds that the ultraviolet divergencies can be absorbed in the coupling constants Ki 
with the following renormalization ot the low-energy coupling constants ||12||: 



A 



K r M 
/i d - 4 



S,A 
1 

16tt 2 U^4 



1 r 



ln(47r)+r'(l) + l 



(3.15) 
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and 



Si — | £2 — Z £3 — — I 

£4 = 2Z £5 = — I £ 6 = |Z 



S7 — S§ — Z £9 — — 1 

S10 = §Z + j En = I £12 — I (3.16) 

E13 = Ei 4 = E15 = 20Z 2 + 3Z+| 

E 16 = -4Z 2 - |Z - 3 E 17 = 2Z 2 -|Z+| 



2' 

C 



Here, the K^(n) are the renormalized couplings at the scale \i. 
The scaling behaviour of K\{p) is determined by the requirement 

^ = (3.17) 

which implies for K\{\x) in the limit d — 4 

Ej , /i 
16^ /io 



%)=W-^M-). (3-18) 



4 Electromagnetic contributions to 77 — > 3n 

The electromagnetic contributions to the decay 77 — > 37r have been considered long 
ago [§, [131 . The transition amplitude is suppressed due to a soft-pion theorem Jl[]. 
At next-to-leading order in a low energy expansion, there is a further suppression 
factor m^/mj^, if linear dependence on the energy of the odd pion is assumed. This 
led to the conclusion that the electromagnetic interactions alone fail completely to 
explain the observed rate. 

Here the problem is reinvestigated in the framework of effective chiral lagrangians, 
as reviewed in section 2. We calculate the corrections of order e 2 p 2 , where p 2 is a 
generic low energy momentum. We work in the isospin limit, as corrections of order 
e 2 (rrid — m u ) are very small. 

As usual, we define the decay amplitude A by 

< 7r°7r + 7T~out I r/> = i (2tt 4 ) 5 4 (P f - P t ) A(s , t, u) , (4.19) 
with the Mandelstam variables 

s = (p v -Pno) 2 t = {jp v - p n +) 2 u = (p v -p^-) 2 . (4.20) 
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Due to G-parity, the three pions emerge in a I = 1 state. Hence, the decay 
amplitude into three neutral pions, A(s,t,u), is also determined by A(s, t, u) 



< 7r 7r°7r°out | r] > = % (2vr 4 ) 5\P f - P { ) A{s, t, u) 

A(s,t,u) = A(s,t,u)+A(t,u,s)+A(u,s,t). (4.21) 

Charge conjugation invariance requires A(s, t, u) to be symmetric with respect to 
interchanges of t and u. 

The electromagnetic contributions to A(s, t, u) may be written as 

Aqed(s, t, u) = e 2 f(A, P 2 , m u , m d , ,m s ,.. .) (4.22) 

where A is the scale of the theory and p 2 stands for any of the variables s,t,u. Keeping 
ratios p 2 :m u :md:m s fixed, / can be expanded according to 

/ = Jo + fi + ■ ■ ■ ■ (4.23) 

fo is of order p°, fi of order p 2 etc. The lowest order term fo vanishes, as was shown 
by Sutherland Jl| . As noted, fi is also suppressed, if the energy dependence of the 
amplitude is linear. This follows from the soft-pion theorem Q 

HmA+i° = limi™ = 0(e 2 m). (4.24) 
Hence, for a linear amplitude 

A£ed = a + b-s, (4.25) 



it follows that at the soft-pion point s = M 2 + 3M 2 , and t 



u 



^i D (Po = 0) = a + b(M 2 + 3M 2 ) (4.26) 
4?ed(Po = 0) = 3a + b(M 2 + 3M 2 ), (4.27) 



and therefore 



m 



0, b = (D( — ). (4.28) 



In chiral perturbation theory the term fi can be calculated in a systematic manner. 
The assumption (f4.25|) is seen to be violated by non-local terms arising from one-loop 
graphs. However, only Kaons appear in the loop, therefore the resulting amplitude 
is smooth throughout the physical region. Equation (|4.25| ) is then still a good 
approximation and the — suppression is effectively at work. This is the main 
reason why the next-to- leading order term fi is small. 

The relevant diagrams are shown in figure 1. The graphs a) and b) are the non- 
local unitary and the tadpole graphs respectively, and c) contains the direct and 
pole graphs. 
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a) 







b) 



< 



c) 



Figure 1: Graphs contributing to the decay rj — > 3n. The filled circles denote vertices 
of order 0(e 2 ) and the filled boxes those of order 0(e 2 p 2 ). 



The electromagnetic vertex in these diagrams is given by the leading order term 
in eq. ( |2.3|) and by eq. ( |3.11|) . There are no graphs with pions in the loop as they 
are in addition suppressed by the 7T° — rj mixing angle e, which is of the order 

O(e)^10- 2 = o C U ~ md )- 
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Also, the mass difference {M\ ± — M^ )qcd contributes only to the order 0( mu ~ md ) 
and can be neglected. 

We then obtain for the unitary corrections 

U(s, t, u) = i2 ^ 6 [3s(3s - 4M 2 K )±C kk (s) + 6(3s - AM 2 K )J kk (s) 

-3(3* - AM 2 K )J kk {t) - 3 (3u - AM 2 K )J kk {u) 

-(5s - AM 2 K ) J kk (0)], (4.29) 



for the tadpoles 



rp/ \ C 9 2Ml 



v/3F 6LV ^4 M 2 -M 2 ' 



2M 2 5s M 2 . . , „ , 

+^-^ + ^]^(0), (4.30) 



and finally for the sum of the pole and direct graphs 
D{s) = -^[1 + 3(S - S0) - 



9V3F 2L Ml -Ml 1 
{-3^ + ^ 4 + K 5 + K 6 -K 9 - K w }. (4.31) 



We have introduced the following functions C kk and J, 



A'/,- 



I Cu(s) = i i ^ arcta , 

i 167T 2 M 2 . a/A VA 



Jkk(s) = — [2 — 2 arctan — =1 + J kk (0) 

K 1 16vr 2[ s JA K ' 



with 



J fcfe (0) = -2A - 2k + 0(d-A) (4.32) 



A = s(4 - s) 



3s = s + t + u = m 2 + 3m 2 , (4.33) 



and 



A = £<dV>< 4 *> + r '« + 1 » 

1 A/f 2 

With equation (|3.15|) and the following definition: 

Jkk(s) = J kk (s) - Jfcfc(O) (4.35) 



9 



we obtain the finite result 

A QED (s,t,u) = e 2 (U(s,t,u) + T(s) + D(s)) 



Ce 1 1 1 I 

^T 6 [^s - AM 2 K )-C kk {s) + -(3s - 4M|) J kk (s) 

~{3t - AM 2 K )J kk {t) - J(3u - AM 2 K )J kk {u)\ 
i ±Ce 2 Ml 3(s-s ) 



with 



+ K r M-K r M-K r M. (4.37) 
Some remarks concerning the structure of this result are in order: 

1. The scale dependence of k is cancelled by the running of K r (fi). Using equa- 
tions ( |3.16| ) and ( |3.18] ) one gets 

3C a 2 

frfr) = K r M - -E^r- 2 H~ 2 )- (4.38) 

2. As expected, the amplitude vanishes in the soft pion point 

s^M 2 , t,u->0 (4.39) 

if the lowest order relation 3m 2 = <Am 2 K is used. This is in accord with Suther- 
land's theorem which predicts a vanishing amplitude if the pion mass is sent 
to zero. 

3. The final result ( 4.36|) consists of a nonlocal piece (first two lines) and a poly- 
nomial to first order in s (last line). The polynomial part is in proportion to 
the current algebra amplitude of strong interactions. It exhibits the M 2 sup- 
pression as implied by Sutherlands theorem. Since the counterterm K r enters 
only this part of the amplitude, it's contribution is suppressed too. The nonlo- 
cal piece circumvents this suppression as it is clearly not linear in s. However, 
since only Kaon loops contribute, this part is kinematically suppressed. 

4. The rate is to a large extent fixed by the amplitude at the center of the Dalitz 
plot. We therefore expand ( [4.36| ) according to 

^4qed(s, t, u) = a + ai(s - s ) + a 2 {s - s ) 2 + a 3 (t - u) 2 . (4.40) 
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Explicitly, the constant term is given by 

«o = |^o(3,o - 4Mi)ick( So ) + 4M* [k - ^K r (^ } . (4.41) 

Note that at s = sq the remaining nonlocal piece in proportion to Ckk is 
numerically less suppressed than for instance Jkk- 

5. The linear slope a\ is given in terms of the same counterterm K r as is the 
amplitude at the center of the physical region. Thus, if isospin breaking in 
the quark masses is switched off, the electromagnetic rate would be fixed in 
terms of the linear slope. The quadratic coefficients 02, get contributions 
only from the nonlocal part of Aq^d- They are given in terms of the known 
coupling constant C . 



5 Numerical Results 

As mentioned, the one-loop result of Gasser and Leutwyler yields a rate of 167 eV, 
much below the experimental value. Although the unitary corrections are avail- 
able H, we will use the one-loop result as reference. This is inessential, since the 
electromagnetic corrections are very small. 

We set F equal to the pion decay constant F n , whose experimental value is F n = 
92.4 MeV if electromagnetic corrections are included [pj| . The low energy constant 
C can be determined by using the underlying vector and axial-vector resonances. 
Using the Weinberg sum rules |f[£l| , one obtains [[L7j 



where M p is the mass of the p-meson, M p = 770 MeV and F p is its decay constant, 
F p = 154 MeV. Numerically, 

§i = 0-84. (5.43) 

7T 

The values of the coupling constants K\ are so far unknown. The 'naive' chiral 
counting law []BJ gives 

I K * 1= 7^ ~ 6 ' 3 x 10 ~ 3 ' ( 5 - 44 ) 

l07T z 

while an estimate for Kl yields [0 

Ktfji, = O.BGeV) = (-1.0 ± 1.7) x 10~ 3 , (5.45) 

consistent with the naive rule ( J5.44j) . 

We conservatively use eq. ( 5.44|) and vary the absolute value of K r between the 
bounds 

I K- |< ^. (5.46) 
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To get a feeling of the numerical size of the elecromagnetic corretions, we expand 
the amplitude (f4.36|) around the center of the Dalitz plot and compare it to the 
ChPT one-loop result for the strong interactions. Writing the Taylor expansion of 
the sum of these amplitudes as 

^4qcd(s, t, u) + A QED (s, t, u) = a + a^s - s ) + a 2 (s - s ) 2 + a 3 (t - u) 2 (5.47) 

we find 

QCD . n.l , I 

= -(0.17 + 0.02z) + 0.14 x 10~ 2 ± 0.30 x 10~ 2 , (5.48) 

where a® CD denotes the strong one-loop value, Oq' 1 the contribution from the non- 
local terms in eq.( |4.36| ), and a l that of the local terms which depend on the low 
energy coupling constants Ki. The variation of a l are due to the changing value of 
K r . Similarly, we obtain 

Ql = a f CD + a? + a[ 

= -(2.25 + 0.75i) + 0.62 x 10~ 3 ± 0.30 x 10" 1 [GeV~ 2 ]. (5.49) 

a 2 = a% CD + a n 2 l 

= (1.33 -0.9H) -0.25 x 10" 1 [GeV -4 ]. (5.50) 

a 3 = ag CD + a$ A 

= -(0.93 -0.41i) -0.85 x 10~ 2 [GeV" 4 ]. (5.51) 



Here, we used everywhere the 'new' value of F n pT| ; the other relevant factors for 
the QCD amplitude are taken from ref 0. From this we see that for the extreme 
values of K r the local contribution dominates over the non-local contribution. Both 
are, however, small compared to the ChPT one-loop result. We therefore expect 
no significant changes from QED to both, the decay rates as well as the Dalitz plot 
distributions for the decays rj — > 7r°7r + 7r~ and rj — > 3n . 

The decay rate r r) _ >7r o 7r + 7r - is obtained by evaluating the phase space integral [] 

dr v ^ n o n+n - = —j^-K p i - P f ) | A QCD (s, t, u) + A QED (s, t, u) | 2 dfi, 

d 3 p 7V o d 3 p 7T + d 3 p w - 
^ (2nY2pl (2n)*2fi + (2tt) 4 2 P °_ " ^ } 

We get for the integrated rate 

r t7 _ >7r o 7r + 7r - = 165 ± 5eV (5.53) 
where the uncertainty is due to our lack of knowledge on K r , and for the ratio r 

1.42 < r = ^ 37r ° < 1.43 . (5.54) 

r\— -+7r 7r+7r~ 
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Note that the phase space is very small and thus very sensitive to the mass difference M 7T ±—M n o 
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In both cases elecromagnetism gives, as expected, small corrections to the results 
obtained by 0, where the purely strong effects have been considered 

r?2S^ = 167eV ( 5 - 55 ) 

r QCD = _ (ggg) 

Even though QED can give rise to a corretion of up to 4 per cent to the decay rate 
T v ^ 7T o 7T + 7T - , the Dalitz plot distribution remains unchanged. This is due to the fact 
that the non-local contributions are small and that the functional dependence of 
the local term in eq.( |4.36|) on the Mandelstamm variable s and that of the current 



algebra amplitude are the same. Since the Dalitz plot distribution is a measure of 
the s-dependence of the amplitude A(s,t, u), the K r (fi ) contribution only leads to 
a different normalization of the distribution, but there is no significant change of 
the slopes. 



6 Conclusions 

In this paper we have calculated the direct electromagnetic corrections to the decay 
?7 — > 37r. Using the general framework of chiral perturbation theory, we determine 
the first nontrivial term of order e 2 p 2 . This includes mesonic loops and a combi- 
nation of counterterms introduced earlier HT2|. Unlike for meson mass differences, 



the corrections are tiny, and amount to at most two percent of the ChPT one-loop 
amplitude at the center of the decay region. f\ The Dalitz plot parameters are also 
rigid with respect to electromagnetic corrections; from the results in section 5 we see 
that both linear and quadratic slope parameters are practically unaffected by the 
electromagnetic corrections. This is clear from the form of the amplitude: First, the 
nonlocal part which circumvents the Bell-Sutherland suppression arises from Kaon 
loops, which are kinematically stongly suppressed. The larger (by a factor of 13) 
pion loops multiply the small quantity (m^ — m u ) and can be neglected. Further- 
more, the polynomial part is suppressed by the small factor M 2 jM\. This includes 
in particular the contributions from the coupling constants JQ. Thus, the depen- 
dence on these little known constants is weak; varying them between the rather 
large values ± results only in a variation of less than two percent relative to the 
one-loop ChPT amplitude of the strong interactions. 

Thus, our result confirms the usual picture, namely that the large decay rate 
requires rather large isospin violations in the strong interactions. Moreover, any 
precise measurement of the rj — > 3tt Dalitz plot yields important information on the 
strong interactions without contamination from the electromagnetic ones. 

Acknowledgements. 

5 There are other, indirect, electromagnetic corrections which influence the fundamental param- 
eters, such as F x - They are known |ll|| and change the normalization of the QCD amplitude 
noticeably, resulting in a 8 % increase of the rate. 
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